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SEPARATOR PLUGS FOR LIQUID HELIUM 


J. M. Lee, S. W. K. Yuan, W. A. Hepler, T. H. K. Frederkmg 

I. INTRODUCTION 


The present report covers the Summer 1984 research. During conferences 
in 1984 it has become clear that previously reported results for porous plugs 
serving as phase separators may be altered significantly by attachment of 
upstream and downstream devices. An example of the latter is a fluid con- 
striction using a heater. The effort of this Summer 84 has been formulated 
originally as a set of three tasks: 

1 . Data analysis 

2. Plug measurements 

3. Critical assessment of applicability of results 

The present report comprises a description of two-f luid-model-relatea permea- 
bility of He II and phase separation in the linear regime (Section II). A 
brief outline of experiments is given in Section III. Further, a critical 
discussion is contained in Section IV. Supplement I presents a conference 
paper dealing with the question of obtaining the Darcy permeability with 
emphasis on "high" temperature methods. Supplement II is a conference paper 
on a related topic, i.e. porous plug use for pumping. These supplements are 
included to give a perspective of "negative" pressure differences, seen from 
the phase-separator point of view. 

II. TWO-COMPONENT EQUATIONS OF THE TWO-FLUID MODEL 

The concept of permeability has been a topic of discussion m various 
fields, such as filtration, chemical engineering-related processes in general, 
vacuum technology, oilfield reservoirs, bio-systems, and others. In recent 
time, the advent of advanced space cryogenics support systems for various 
tasks has added new phenomena at low temperatures (e.g. 1.8 X) . The use of 
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porous media, such as sintered metals, with nominal pore sizes extending from 
the order of 0.1 to 10 inn has posed several questions. Major transport 
mechanisms have been clarified in the past few years, however much remains to 
be done m order to get a comprehensive picture. 

The present work continues preceding investigations in the area of space 
cryogenics related to the task of venting and separating vapor from liquid 
Helium II. The generally non-Newtonian transport is expressed in terms of the 
Darcy permeability for Newtonian fluids in the Stokes-Darcy regime. Recent 
results confirm that "permeability" is indeed useful when the unique defini- 
tion of the Darcy permeability is kept m mind. Uniqueness requires the lim- 
iting value, as the speed goes to zero, and exclusion of Knudsen transport. 
Further, opening and closing of "dead space" pores is excluded as well. 

The emphasis is on the evaluation of various effects which modify the 
permeability. The two-fluid equations are used. The interdependence of vari- 
ous rates of transport during phase separation is considered. One-dimensional 
equations are used. 

The linear regime of transport is discussed for a vapor-liquid phase 
separator (VLPS) based on the use of the fountain effect in He II. The 
thermo-mechanical pressure is directed toward the interior of liquid He 1 *. The 
liquid vessel is below the lambda temperature of the He I-He II transition. A 
low temperature T is on the vapor side and T+AT is on the liquid vessel side. 
Thus, liquid is kept inside the tank when design conditions are met. Various 
assumptions underlying a basic system of equations are discussed with subse- 
quent formulations of vector and scalar quantities. The results are inspected 
m the light of findings reported in the literature, and m the light of 
experimental evidence. 


Device and System Constraints 


Consider a phase separator during steady transport of entropy, liquid 
mass, normal fluid, superfluid, and heat. The five basic statements are the 
density postulate for the two-fluid system, the mass flux density postulate, 
the heat flux density, the relative velocity, and the first law constraint for 
the vapor-liquid interfacial domain at the downstream end of the separator 
(Fig. 1). Downstream and upstream is viewed from the perspective of mass 
throughput and heat flow. 

Basic two-fluid density postulate . The density of liquid helium is 

P = P 3 + Pn (1) 

(subscripts n and s denote normal and superfluid respectively). The usual 
two-fluid approach ignores the condensate fraction referred to in approaches 
aiming at contact with quantum mechanical parameters (e.g. Ref. 1). 

Basic two - fluid flow postulate . After the introduction of postulate (1) 
it appears to be plausible to assign velocity fields to the two fluids assum- 
ing that, at low energy, there exists independent motion to first order: the 
mass flux density "j is the sum of two contributions: 

^ ^ -J*. 

J = pv = P s v s + p n v n (2) 

In the formulation of the related quantum hydrodynamics it has been considered 
useful to introduce a new frame of reference moving with the superfluid velo- 
city v This frame replaces the Euler frame adopted in most cases (e.g. 

Ref. 2). For instance, the mass flux density (momentum j) is written as j = 

pv s + Jo ( J o 1S the J" va l ue in the new frame "K"). The velocity v“ of a given 
liquid element is zero m the frame "K" , and the liquid moves with the rela- 

-i. 

tive velocity (v^ - V g) in "K". 

Relative velocity w. In the preceding remarks, the Galilean principle 
of relative motion has been referred to. Because of the general utility of 
this type of approach, we employ the relative velocity subsequently: 
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LIQUID 



Figure 1. Schematic diagram of porous plug operation 
in He II in the phase separation regime ; 

V g Superfluid velocity ; 

V n Normal fluid velocity ; 

V Mass flow velocity . 




( 3 ) 


Heat flux density (q) . Original versions of the two-fluid model have 
been proposed by London and Tisza, and by Landau. Usually, the starting point 
has been the assignment of entropy to the normal fluid. Thus entropy convec- 
tion and heat convection is brought about by flow of normal fluid. This leads 
to the heat convection statement (e.g. Ref. 3)- In addition, London 1 * formu- 
lates (pS) = (P n S n ) with S = entropy per unit mass. At very low temperature, 

T < 0.6 K, phonons are the dominant entropy carriers. The thermodynamics of 
ideal phonons leads to the results that (p n /p)/(ST) = const = Op/9P) s (i.e. 
reciprocal square of the first sound speed); (P pressure). Thus, we have a 

limiting value of iim s = liml— | = =lim JL 00 The normal fluid entropy 

T-o n t-kAV T-0 ST 

appears to grow beyond all limits as 0 K is approached. No further use of 
this relation will be considered subsequently. Instead, in the frame "K", 
mentioned above, the heat flux density q is expressed as 

q = P s wST (4) 

First law constraints at the vapor - liquid interface . The first law of 
thermodynamics for steady flow has been considered frequently. According to 
the simplest version, the net change of the enthalpy of the flowing streams, 
across a control volume, bounded by suitable control surfaces, is the sum of 
the shaft work supplied and the heat transferred. At the ideal interface, no 
shaft work is made available. Further, per unit area, we have a heat input of 
q on the liquid side of the interfacial domain (e.g. Ref. 5). Close to the 
interface the liquid has a low order parameter. Therefore the classical ver- 
sion of the first law is expected to be meaningful. The enthalpy flux enter- 
ing is m^n^ while the flux leaving is m y H v . Thus, the net difference is 
readily formulated. For steady flow, we have continuity m mass transport, 
i.e. m = m = m. Therefore, we abbreviate j = mass flux density (m/A); 

Lt » 

j(h v -h l ) = j* = q (5) 
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The enthalpy difference (H y -H L ) is the latent heat* of evaporation (X). Con- 
cerning further details of interfacial studies we mention a recent investiga- 
tion of low pressure evaporation (G. Bertrand and R. Prud ’Homme^) . In addi- 
tion, H. Wiec'nert^ has studied details of the boundary conditions for the 

O 

liquid-vapor interface of He II. Schotte and Denner introduce kinetics 
effects arising when the mean free path becomes important, i.e. for the low 
pressure domain For further information, these references may be consulted. 

in summary, the present modeling approach has five basic equations which 
are useful for a better understanding of the transport processes at least in 
the linear regime of VLPS systems. 

Vector quantities . Consider the five basic equations listed in the 
preceding discussion. The normal fluid may be eliminated from Eq. (2) using 
Eq. (3): 

-w -a* -a* 

J = v 3 p + wP n = q/* (6) 

Elimination of w from Eq. (6) by means of Eq. (4) results in 

q = -v s p s ST/[(p n /pJ - (p s ST/Ap)] (7) 

In the vicinity of the lambda temperature (T.) we have p /p -*■ 1; (p /p) << 1. 

a ns 

The expression in brackets on the right hand side of Eq. (7) is positive. It 
turns out that this is also the case at all T in the roton-dominated range. 

J*. 

Thus, the heat flux density q is convected m counterflow to the superfluid 
motion . 

The expression [(p / p ) - ( p / p ) (ST/ \ ) ] is shown in Figure 2 
versus T. 


* H y — h, = X only if AT is zero. At the interface the temperature of the vapor 
and liquid are the same for infinite radius of curvature and high pressures of 
continuum regimes (we refer also to Appendix C). 
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Figure 2. Flow rate factor as a function of the 
temperature; Equation (7): q/(-v s p s ST). 
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Similarly, the superfluid velocity may be eliminated from Eq. (2) using 
Eq. (3)* The equation for j is written now as 


J = pv n ~ P s w (8) 

Again, the relative velocity is eliminated using Eq. (4): 

-k. 

q = + P v n STA/( A+ST ) (9) 

The expression in the denominator of the right hand side of Eq. (9) is posi- 

-fa. 

tive. Thus, q is in the direction of v nf in agreement with normal fluid 
behavior postulated for the two-fluid system. According to eq. (9) one may 
write 

q = q ZNMF (1 + ST/A) -1 (10) 

The reference quantity is the zero net mass flow rate q ZNMF . As T is lowered 
sufficiently, ST becomes much smaller than the latent heat. Therefore, the 
solution for zero net mass flow is approached asymptotically, as discussed 
elsewhere^. Details of the preceding treatment are in Appendix A. 


Mean transport rates and mean speeds . The special flow constraints of 
the zero net mass flow (ZNMF) mode are rather well known. There appears to be 
no rigorous local condition derived from Eq. (2) for j = 0. In other words, 
the flow fields of Newtonian-like normal fluid flow with a parabolic velocity 
distribution, and those of the superflow with a "plug flow" distribution per- 
mit satisfaction of j = 0 only at one particular location of pipe flow. How- 
ever, mean values of the flow equation do satisfy the scalar equation 

l J sl = Mnl’ or 

l J sl = P s I v s I = p n I v n I = Uni (ZNMF) (11) 


Figure 3 represents examples of vectors compatible with Eq. (2) (case 
(a) is a possibility which might be realized m a porous plug; case (b) is the 
oarallel vector arrangement). The inset of Figure 2 presents the parallel 
case for the relative velocity (w), the superfluid velocity and the mass flux 
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Figure 3. Examples of vectors of the two-fluid model description 
of vapor-liquid phase separation; 

a. Different directions of superfluid velocity .normal 
fluid velocity and mass flow velocity ; 

b. Parallel velocities . 
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density j . 


The mean heat flux density of Eq. (10) is written as 

M = +p s I v s I ST/[ (p n /p) ~ (p s /p)(ST/l)] (12) 

Some additional details are given in Appendix B. 

Van Sciver 10 has presented an interpretation of the transport phenomena 
resulting from a superposition of a transport velocity on the specific case of 
Gorter-Mellmk counterflow. The underlying equation is deduced from an adop- 
tion of an Euler coordinate frame. Thus, the usual thermal energy equation is 
replaced by 

p DH/Dt = p[3H/3t + v3H/3x] = -3q QM /3x (13) 

“ q GML * -f (T)( 3T/3x) 1/3 ] is the Gorter-Mellmk equation for the counterflow 
heat flux density; f(T) is the Gorter-Mellmk temperature functional; H is 
enthalpy. For steady flow we have 3H/3t = 0. The preceding one-dimensional 
equation for flow along the x-direction becomes 

dH/dx=+f(T)d(dT/dx) 1/3 /dx ; AT << T (14) 

The previous transport equation of the two-fluid approach is not entirely con- 
sistent with the special equation (14). Once consistency with Van Sciver's 
approach is incorporated, the preceding equation for q (Eq.10) is to be 
replaced by 

q ’ "ZNMF* 1 - STArl <l5) 

The factor (1-ST/A) -1 is displayed m Figure 4 versus the temperature. The 
two-fluid result, Eq. (10), is included. The difference between the two-fluid 
mode Eq. (10) and Eq. (15) is about 23$ at 2.1 K. 
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HEAT FLOW RATIO q/q 



Figure 4. Heat flow ratio of ZNMF as a f unct ion of the 
temperature . 
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Discussion of T- influence and related work . In general the superposi- 
tion of mass flow and heat flow may be arranged Independently by external 
means. For the special geometry of transport through a duct with adiabatic 

^ -Jk 

walls one may distinguish the two cases of q and j in parallel, and q and j m 
counterflow, in the VLPS system however, the first law constraint (Eq. 5) cou- 

— »» Jk 

pies the two transport rates of heat and mass. Thus, q and j are kept in 
parallel. The analysis of Van Sciver 1 ® has been based on the turbulent flow 
studies of Johnson and Jones. For the parallel case there is an enhancement 
of q beyond the ZNMF value. In contrast, the two-fluid model Eq. (10) 
predicts a lowering of q with respect to the ZNMF mode. 

It is noted that the idealized vapor-liquid interface may have a less 
ideal domain than presumed m the linear system of equations. This is 
expected when turbulence and/or vibrations start to shake the interface. How- 
ever, "classical" findings on evaporation in the presence and absence of mono- 
layers have resulted m the conclusion that pure heat conduction in the 
interface-adjacent layer is the dominant mechanism (G. T. Barnes and D. S. 
Hunter 11 ). Thus, for a certain class of VLPS phenomena, a similar restriction 
related to He I-like diffusion is expected to dominate the liquid layer of low 
order parameter close to the vapor. Concerning details of the microscopies, 
we refer also to the review of R. Evans 1 ^ on "The Nature of the Liquid-Vapor 
Interface". The He II-vapor interface has been treated by Hills et al. in 
another paper 1 3. 

The temperature difference due to surface tension is outlined m Appen- 
dix C. The one-dimensional analysis of the non-lmear transport regime by van 
Sciver is sketched in Figure 5 versus the flow parameter used in Ref. 10. 

When the entropy and mass flow rates are in parallel as for the VLPS mode, the 
heat transport rate is increased. It is noted that this increase m the aata 
range documented by Johnson and Jones 111 does not exceed a factor of 2. As 
discussed previously, the vapor-liquid phase separator has entropy flow m 
parallel with mass flow. 
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Figure 5. Heat flux ratio vs. flow parameter 

K ys of Van Sciver (ref. 10); 

K vs is proportional to the thermal 
energy convection rate of forced 
flow with speed v to the Gorter- 
Mellink convection rate. 



Another example of indirect evidence supporting a lowering of the heat 


transport rate during counterflow of entropy and mass has been presented by 

1 C 

Tough and Baehr . Critical conditions have been considered. Further, it is 
noted that attention is on the duct downstream of a pump and not on the porous 
medium. A fountain effect pump is used. It causes flow of liquids in coun- 
terflow with the normal fluid flow in the porous medium. At critical condi- 
tions, characterized by the subscript c, upon an increase of | v nc j, the flow 
rate expressed as |v | i s decreased. Thus, there is partial confirmation of 
the consequences of the Van Sciver approach for the critical conditions of the 
two-fluid system. The results of Tough and Baehr 1 ^ are given for 1.4 K and 
1.6 K. The data are contained m Figure 6. 

III. EXPERIMENTS 

The apparatus used has been similar to 
The change incorporated in the present setup 
flow modulator studied previously. An outer 
central vent tube with a porous plug located 
tube. Vent conditions are modified by a set 
terfly valve ) . 

Procedures are identical as far as the general preparation is concerned. 
During the tests, various positions of the valves have been used to cause dif- 
ferent vapor flow impedances. The range of flow is varied in a finite mass 
flow rate interval caused by heat supply to a heater located in the "bath". 

The bath proper is a jacket surrounding the central tube. The jacket in turn 
is immersed in the outer bath. The liquid level in the jacket is a monotoni- 
cally decreasing function of time during the tests caused by the heat supply. 
After several runs the liquid level oecomes too low. Once tms state has been 
reached an outer "thermo-pump" is activated by heat supply at the downstream 
end of the pump. The liquid is transferred into the jacket by this pump using 
the fountain effect, as described in Supplement II. 


the one described previously 1 ^, 
has been the elimination of the 
liquid Helium II bath surrounds a 
at the lower end of the vent 
of two valves (bypass and but- 
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Figure 6. Results for counterflow of mass and 
entropy in a porous medium for 
critical conditions (Tough-Baehr ^ ) . 


Note: Tough and Baehr place emphasis 
on a duct attached to the pump 
using negative flow velocity as 
sign convention for this duct. 
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The resulting flow rates deduced from level changes are affected by 
gross unbalances m heat supply, bath temperature, and valve position. By 
varying valve cross sections available, in order to match the heat supply 
rate, a certain "bath" temperature interval is covered rather close to quasi- 
steady conditions. (If this is not the case the bath temperature either 
drifts down or increases upwards).. Data taken are affected by the degree of 
stability created by the valve impedance selected. The ideal case of constant 
mass flow rate is enforced when there is a large impedance. However this case 
is characterized by rather small throughput of mass through the plug. There- 
fore, a number of runs are established at conditions farther away from the 
constant mass flow case. This is necessary in order to cover a reasonable 
mass flow range in one set of experiments. 

Figure 7 shows Darcy permeability results for various plugs at room tem- 
perature. Dg is the equivalent Ergun diameter. This diameter has been modi- 
fied using the porosity e. The throughput length L is the square root of the 
Darcy permeability. If there exists a packing geometry exactly equivalent to 
a spherical particle system rather closely packed, the diameter D £ is expected 
to vary but little. The actual data spread indicates a finite spectral range 
associated with the particular plug system under consideration. Details of 
the underlying concepts are given m Supplement I. 

Experimentally obtained quantities are displayed m the next Figures 8 
through 13. The variables are shown versus the time t as obtained during the 
experimental runs. 

Figure 8 presents the decrease of the liquid level z with time t. From 
the conservation of mass condition, velocities and thermal energy convection 
rates through the porous plug are calculated. The mean velocity is displayed 
in Figure 9. The heat flux Q is obtained from the first law constraint at the 
plug exit (Figure 10). 
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Figure 7 . Distribution of the porous plugs classified 
by the Darcy characteristic length. 
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Figure 8 . Liquid level change with time during VLPS run; 

(Plug : Bronze; nominal pore size 5-15 jum) ; 
(This plug has been used in the various runs 
for which data are shown in subsequent figures). 
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The bath temperature change during slow changes of state is plotted as a 
function of time in Figure 11. The bath temperature is the upstream tempera- 
ture T u of the plug. Figure 12 shows the temperature difference across the 
plug. Figure 13 presents the pressure difference obtained from the pressure 
transducer readings. 

Once steady transport or quasi-steady conditions are secured, the 
"effective permeability" is readily available. In general, the data are in 
the non-linear regime for the plug results of the preceding Figures 8 through 
13. 

The mass flow rate m expressed as v q 1S shown in Figure 14 for three 
different bath temperatures of 1.55 K, 1.6 K and 1.85 K. 
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Bath temperature as a function of time 








IV. DISCUSSION AND CONCLUSIONS 


The present set of equations may be used to find the mass flow rate m 
versus the pressure difference AP y in the linear regime of the porous plug 
phase separator. In the analysis only a few equations are needed to obtain 
the mass flow rate from an integration. The following equations are used for 
steady transport at relatively large AT and AP y respectively: 

1. Coupling condition for entropy and mass flow. 

2. Mean normal fluid convection rate as a function of the driving tempera- 
ture gradient. 

3. Ratio of mean q Q -value of the porous medium to the mean q Q of zero net 
mass flow. 

4. Mass flow rate as integral of a T-dependent function F(T) (presented 
below) . 

5. Vapor pressure function (T-58) for conversion from AT to AP y . 

The preceding set leads to the mass flow rate of steady transport m terms of 
thermophysical properties: 

T » 

a ' K pn<*tot /L > { F<T,<1T (16) 

K pn is the normal fluid permeability (Darcy law analog in the linear regime). 

A is the total plug cross section. L is the axial length, i.e. the depth 
of the plug. The property integral is written as follows: 

F(T) = p 2 S 2 T(1-ST/A) -1 n" 1 A -1 . (17) 

(p liquid density, A latent heat of vaporization, n n shear viscosity of the 
normal fluid). 
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Figure 15 shows the ratio of (m/K^) as a function of temperature when 
the T-range covered extends from a bath temperature of T all the way to the 
lambda temperature T^. For any other interval (T 2 ~T-|), the mass flow rate may 
be determined using the difference function 



( 18 ) 


Figure 16 presents (m/K^) as a function of the vapor pressure differ- 
ence AP v . The function is qualitatively similar to Equation (16). In both 
Figures 15 and 16 a ratio (A^^/l) of unit length has been introduced for sim- 
plification. 

Use of the thermophysical property integral for an arbitrary transport 
case usually leads to a non-linear Kp n -function. In experiments, however, 
this may also indicate that the linear regime was not established over the 
entire length of the plug. A VLPS example is shown in Figure 17. The K pn - 
value obtained at the lowest temperature is much closer to the Darcy analog 
value than at any other temperature. The data constitute "effective normal 
fluid permeabilities". It is noted that the (dashed) reference value is for a 
similar plug (1 pm nominal pore size, material: aluminum silicate). (Further 
details are given in Supplement I.) The phase separation process is quite 
close to the ZNMF mode. Therefore another comparison is useful for under- 
standing of porous media observations. 

Figure 18 is a data set for the ZNMF mode of the early pioneering work 

of Allen and Reekie 1 ^. This particular result has been reproduced in the book 

1 ft 

by W. H. Keesom, Liquid Helium . A powder-filled tube had been usea by Allen 
and Reekie. Therefore, this set of data appears to be among the earliest 
porous media results for He II. The original figure of these authors (Fig. 2 
of Ref. 17) appeared to have been a matter of considerable discussion 
because of the use of the "total conductivity" of the packed bed, given as 
"(watts/K.sec)". Professor Allen kindly responded to an inquiry about this 
point with letter of June 25, 198^. He notes that "it is hard to recall the 
precise way m which the measurements were made m the paper with the late Jim 
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Figure 15. Normalized flow rate to permeability- 

ratio for (A^ ^/L) of unit length (1 cm) . 

Integral for T-range from T to the lambda 
temperature is plotted as value at the low 
initial temperature on the downstream end 
of the plug. 

Numerical values have been selected m 
order to compare with ref. 20 (ly.m plug). 
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Figure 16 . Normalized mass flow rate to per- 
meability ratio of Fig. 14 vs. 
pressure ; 

The integral all the way from the 
downstream state to the lambda point 
(saturated liquid) is plotted at the 
downstream value of P. 

Numerical values have been selected 
m order to compare to ref. 20 ( lum 
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Figure 17 . Effective permeability K of flow 

pn 

through 1 pm plug (norminal size, of 
Klipping group (20) ;) (Evaluated: Cour- 


tesy: Chau Nguyen): 


Note: the data trends are reminiscent 
of the early results of Allen- 
Reekie (ref. 17) in the non- 
linear regime (Figure 18). 
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Reekie, especially as my old note books do not seem to have survived. How- 
ever, re-reading the description of the measurements of Figure 2, I think the 
best description is in the first three lines at the top of p. 118, which 
refers to Figure 1(b). It was the wattage put into the heater H divided by 
the rise in temperature at the thermometer T when equilibrium was reached. T 
had been calibrated against the vapor pressure, and this was also measured 
when equilibrium was reached. As far as I recall, we did not have a second 
resistance thermometer in the bath, but we did have a good cathetometer . " 

Therefore, the letter indicates rather convincingly that the units might 
have been "watt-sec/ (K sec)", i.e. the thermal conductance [Q/AT] is given as 
[W/K]. The Allen-Reekie result in Figure 18 has been reproduced m this form 
as Q/AT = A tot q o /AT. At high values [Q/AT] at 1.9 K there is a very non- 
linear behavior. At the lowest bath temperature however the data appear to 
produce a nearly constant conductance. In other words, the data are close to 
the limit of K pn = const. 

Conclusions . The following conclusions are suggested from the present 
Summer. 1984 studies: 

1 . The VLPS results appear to be maeed close to zero net mass flow in 
agreement with theoretical prediction. 

2. There is a small enhancement of the heat flux density q Q with respect to 
the zero net mass flow mode predicted by the Van Sciver model. Experi- 
mentally this enhancement seems to be masked by other secondary effects. 

3. The results obtained with several plugs appear to be consistent with 
zero net mass flow in so far as there are considerable departures from 
the linear regime prediction. This is confirmed for the present plug in 
the pore size order range of 10 ym (nominal value). 

4. The normal fluid analog of the Darcy permeability is usually an "effec- 
tive" value outside the linear regime. Procedures to come close to the 
Darcy analog case are most successful for large pores when the 
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temperature is lowered and when AT is decreased. 
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APPENDIX A 


FURTHER DETAILS OF VECTOR QUANTITIES OF THE TWO-FLUID MODEL 


Consider the five equations for the superfluid He II density 
system, for j , for q in terms of the latent heat "X , for q in terms 
of w , and for w : 


? ‘ ?b + Pn 

J-jW. 

= q /> 

q = P s w S T 


w = V - V 
n s 


Elimination of v n to obtain q(v g ). Use of A. 4 and A. 5 

v_ = w + v„ leads to 
n s 

j = ?n w + ? v s 

Elimination of w, employing A. 4 and A. 6 .yields the 
function : 


I 


q = - P 


v S T- 

S U 1 - <Ps/Pn> 


(A.1) 
(A. 2) 

(A. 3) 
(A. 4) 
(A. 5) 

in A. 2 with 


(A. 6) 

desired 


(A. 7) 


Note : The denominator of (A. 7) is always positive in the range of 
roton excitations, and q is opposite to v g . 

Elimination of v to obtain q(v ) . Use of A.l and A. 5 in A. 2 
s ^ n 

with v = v - w leads to 
s n 


? v n ~ P £ 


w 


(A. 8) 


The relative velocity w may be eliminated again using the above pro- 
cedure. The resulting function for q is 

' ' p SI <A ' 9> 


The factor with the latent heat "X is always positive, and q is 
in the direction of v n , as required by the 2 nc * q aw 0 f thermodynamics. 
For illustration, vector examples are shown in Figure A.l 
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(a) Arbitrary case 






Superfluid velocity 
component opposite to 
mass flow component 


(c) Similar to case (b) 

"3,, with small mass flux density 




Figure A.l . Various vectors 


36 


APPENDIX B 


MEAN TRANSPORT RATES OF THE TWO-FLUID MODEL 

Consider mean values of the previous two-fluid equations (e.g. 
Appendix A). For the phase separator processes, the first law and the 
second law of thermodynamics requires that j and j have a flow 
direction from the liquid vessel to the vent lines system. , i.e. 

J j ! + |j s | = | j n J . The preceding equations of Appendix A are re- 

written as follows: 

1 j! + 1 j s ! - 1 j n l (B.D 

1 j J = I qj /X (B.2) 

I q I = p s I w I S T (b.3) 

The preceding equation (A. 6) of Appendix A becomes 


I j I + pK! = Pn I W 


(B.4) 


Elimination of w using (B.3) and (B.4) leads to the 
function [q ( [ I vj]: 


[qj * p |v s | S T 


The mean heat flux density 
|v n | noting that 



is obtained readily 


(B.5) 


as a function of 


I v s I = I w I - I v n J (B.6) 

Further, the mass flux density is 



The resulting heat flux density is 

I q I = p| v n |ST (B.8) 

' " Ov+ s T ) 

Figure B.l is an example of the relative velocity w in 
relation to the other velocities of the two-fluid model. 
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APPENDIX C 


TEMPERATURE DIFFERENCE ACROSS THE VAPOR-LIQUID INTERFACE 

The temperature difference arising from surface tension of 
a curved interface is not zero. Consider an interface characterized 
by a radius of curvature R . The pressure difference associated 
with a spherical shape is 

AP - 2 s' / R (C.l) 

A first order approximation of the related temperature difference 
is obtained from the Clausius-Clapeyron equation for the vapor 
pressure. At the low pressures P of the He II- vapor equilibrium 
the AP V along the vapor pressure curve and AT are related by 

AT = AP v T /( p v > ) ; AT « T (C.2) 

After insertion of (C.l) with AP = AP V into (C.2) we obtain 
a temperature difference of 

AT - ( 2 V /R) [T/(p v > )] : ( AT C<T T) (C.3) 

( surface tension = 0.33 erg/cm 2 at 2 K ). For instance , at 

2 K, a radius of curvature of 5 pm gives rise to a temperature 
difference of approximately 14 milli-K . 
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APPENDIX D : APPLICATION OF POROUS PLUG FLOW 


RESULTS FOR FOUNTAIN EFFECT-BASED 
REFRIGERATION SYSTEM NEAR 2 K 


The system including performance figures for a fountain 

effect refrigerator have been discussed for one particular version , 

called the He II vortex refrigerator , by the originators 

22 ) 

Staas and Severijns . The present work has given considerable 

insight into various operational characteristics associated with 
fountain effect use in porous plugs. A simplified, ideal cycle is 
discussed . 

Consider the ideal cycle based on the sequence of changes of 
state as follows : p - P - p-P ;(P = pressure; p = chemical potential). 

The ideal cycle is shown in Figure D.l in the InT-ln S -diagram , and 
in Figure D.2 in the In P - T diagram. Ideal heat exchange is 
presumed, and ideal "superleaks" . It is noted that not all aspects 
depicted here schematically have been clarified at this time. In 
particular we remark that the fountain effect in conjunction with 
the increase of S with a rise in pressure causes peculiarities in 
cycle operation not known in other refrigeration cycles. Nevertheless, 
local thermodynamic equilibrium is presumed such that one can draw 
continuous curves for the changes of state considered. 

At point 1 a reference state is attained,e . g. saturated 
liquid He II. The change of state from 1 to 2 is a compression 

from to a higher pressure ?2 • The change of state from 2 to 3 is 

assumed to be an ideal isobaric heat exchange path . For instance , 
an aftercooler may operate from to during this change of state. 

If a rather low T is desirable , a counterflow exchanger may be 
employed for the subsequent lowering of T to T^ . The change of state 
from 3 to 4 is ideally at constant chemical potential . This is 
similar to the ideal pressure increase from 1 to 2 , however at 3 
the pressure is being lowered to the initial pressure P^ = P^ . 
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At point 4 the refrigeration load is taken up by the working 
substance up to . If the counterflow heat exchanger is left 

out , the system is modified, as sketched in Figure D.3 . In this 
case the refrigeration load is modified to include the temperature 
range from to . 

It is noted that power cycle operation is restricted in the 
amount of work to be gained from an available heat . The volume 
of the He II changes but little during the cyclic operation. 
Therefore, a rather small amount of work is extracted from each 
cycle. In contrast , the FER operation is characterized by a 
reasonable coefficient of performance once heat exchangers can be 
compact enough with extended surfaces such that a large area is 
available for a relatively small volume , e.g. sintered exchangers. 

An estimate of the performance is given based on a simplified 
result for the net work. It is noted from Figure D.4 that the density 
varies insignificantly with T, and it is primarily a function of the 
pressure. As the cycle takes place between two pressures , a first 

order approximation for the net work is AP-A(l/p) ; p = density. 

2 2 
As we have A ( 1 /p) cs Ap/p > the net work is about AP Ap/p 

— B 

The isothermal compressibility K™ is of the order of magnitude 10 
2 1 

cm /dyne ; (K^ = p (b p/ bP)j, ) . Thus, we have Ap » p AP , 
and the net work is 

<^> dW a ( AP/p) Kj. AP = W (D.l) 

The refrigeration load of the simple cycle D.3 is c" (T^-T^) . 
Consequently, for the simple estimate based on Equation (D.l) we 
obtain a coefficient of performance of 

Q c / w = <VT 4 )/ W (C.2) 

3 _ 

For AP of the order 0.1 bar , with p ~ 0.1 g/cm , for c p of 

the order 1 J/( g K) , (T^-T^) of the order of 2 K , one 

obtains Q /W of the order of 10* ; ( Kapitza resistance 
excluded, W/°0.1 J/kg ; Q ^1 J/g ). 
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COMPARISON OF FLOW THROUGH POROUS PLUGS : 

IS THE DARCY PERMEABILITY A USEFUL MEASURE OF THROUGHPUT ? 

J.M. Lee* W.A. Hepler, S.W.K. Yuan and T . H . K . Frederking 
University of California, Los Angeles , CA 90024 , USA 

ABSTRACT 

Size characterization of convective transport through porous plugs 
in the past appears to have failed to show meaningful lengths gen- 
erally accepted for the description of throughput in superfluid 
He II - related applications. It is the purpose of the present 
studies to obtain the characteristic throughput length L c (=square 
root of the Darcy permeability ) , and to evaluate various phe - 
nomena causing departures from Darcy's law. The nominal size range 
studied extends from the order of magnitude 1 um to 10 um. Data in- 
dicate that L , in conjunction with the porosity , is a quantity 
useful for evaluation of the statistical distribution of manufac- 
turing-related plug spectra. 


INTRODUCTION 

The quantification of flow of cold fluids through porous media has 
received attention recently in conjunction with vapor liquid phase 
separation components for liquid He II vessels operating in space 
/!/ . In the past sintered plugs have been used most frequently 
for filtration tasks. Parameters pertaining to exterior surfaces of 
the plugs have been of prime interest. They may depart from interior 
properties dominating the throughput . The latter ought to be known 
for optimum liquid retention m the vessel, at minimum liquid He 
II losses, by means of fountain effect forces. 

In general, two-phase conditions have to be addressed. However , 
single - phase flow appears to be known only to a limited extent for 
the sintered metal plugs of interest. Recent studies /2/ indicate 
that the length L may be uniquely defined, provided the plugs are 
clean, and the transport takes place in the Stokes -Darcy regime at 
small laminar flow speeds. Several parameters may cause departures 
of the effective permeability (K gff ) from K D . Therefore, the 
present work is concerned with phenomena disturbing Darcy con - 
vection. It is noted that "permeability” in the literature often 
is used as an operational term. An example is the permeability of 
diffusive permeation of gas through glass dewar walls, e.g. Dush- 
man /3/. In contrast, the Darcy-value K Q is discussed in the pres- 
ent work. The experiments are outlined , and finally conclusions 
are drawn concerning the usefulness of . 

CONVECTION PERMEABILITY OF SINGLE-PHASE FLOW 

In sintered porous plugs, the well-defined conditions of monosized, 
near-spherical particle assemblies of packed beds are approached 

^ Present address: NASA AMES Research Center , Moffett Field, Calif. 
+) Presented at Space Cryogenics Workshop, Berlin, Aug. 1983. 
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only in restricted cases. Nevertheless related equations, such as 
Ergun's function, have been used as reference . In addition, the 
simplest geometries may provide ideal permeabilities for infinite- 
ly thin walls in the Stokes regime of low speed. According to 
Darcy’s law_, the volumetric flow rate, per total cross sectional 
plug area, (v G ) is proportional to K D , to the driving pressure 
gradient, (grad P) , and inversely proportional to the shear vis - 
cosity (>j) . 

v Q = K q I grad P I / »j (1) 


For a fluid-filled annulus , with diameters D Q and D 1 , the Darcy 
permeability for very thin walls may be expressed as 


K- 


D eff Z 32 


( 2 ) 


with D 
zero inne 




= <D 2 +D. 2 ) - (D i 2 - o/)Z In (D /D ) 
diameter the tube result (D 0 2/ °32 ) 


2 - 2», , i n ^^ mlt Q f 

iciu iiinci uxoiuc lc; wiie uuut: m&uit \ uq -/ j is recovered , 

and in the limit of a very narrow slit , the value (D 0 2-Dj_2)/ 4 g 
is obtained. 


Among the single-phase phenomena which may perturb the flow are 
the following conditions: 1. Temperature effects; 2. High-speed phe- 
nomena of incompressible fluid flow; 3. Compressible fluid flow ef- 
fects; 4. Mean free path effects during transition toward Knudsen 
transport. 

Temperature influence . In the present experiments with gas, the 
mass flow rate has been measured downstream at room temperature 
T = T e when the fluid has the density p = p . The resulting value 
of K d at low speed and T 4= T g is 

k d = <p e /p> v o<Te H<T)/ {grad p| (3) 

Thus, any T-error caused, for instance by walls at T lower or higher 
than the fluid temperature results in an effective value K 
smaller or larger than K Q . 

High - speed effect . For incompressible fluid, inertia may become 
significant such that vortex-shedding processes are initiated 
during the transition into the turbulent flow regime. Ergun's equa- 
tion Z 1 ^/ incorporates 0. Reynold's linear superposition /5/ of a 
laminar (Stokes-Darcy) term and a v 0 - contribution to grad P . 
The resulting effective permeability is 

K eFf = K D /[ 1+Re 1.75/(150 3 ) 1/2 ] (4) 

(L c = F d 1 Z 2 ; Reynolds number Re = L c p v Q / ; z porosity) . 

Compressible flow effects . There are limiting boundary conditions of 
isothermal and adiabatic flow. An approximation for isothermal flow 
of an ideal gas has been evaluated omitting the Mach number of com- 
pressible flow, e.g. /6/ . For large pressure ratios, (P on the up- 
stream side to P on the downstream side), beyond the critical value 
for a nozzle of about 2, local attainment of the sound speed is 
probable at the plug exit. The isothermal K s ^ c is known to be above 
Eq. (4) . Thus, the adiabatic K e ff is expectea‘to be even larger .As 
little residence time is available at hign speed, the fluid may be 
close to the adiabatic case. 

Knudsen transport influence . When the mean free path starts to be- 
come comparable to the pore size of the plug, particle collisions 
with solids are less frequent than in continuum flow. The resulting 
| grad P| is smaller for a given flow rate. Thus, K e ff becomes 
larger than . 
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EXPERIMENTS 


Experimental runs with He^ gas have been conducted with the setup 
shown schematically in Figure 1. Most of the data have been obtain- 
ed at room temperature with some extensions to tne liquid Nj boil- 
ing point. The gas is admitted from a cylinder, via a pressure re- 
duction valve, to a liquid N 2 trap, transferred into a liquid N 2 - 
dewar and passed through the porous plug. Mass flow rates are 
measured downstream. In the course of the work the apparatus, in- 
cluding the thermometry , has been modified. Fluid temperatures 
are measured upstream and downstream of the plug. Data showing 
significant T-discrepancies are omitted. Cold fluid data are con- 
veniently taken during slow warmup of the system toward room tem- 
perature. Pressure differences are measured with a Validyne P - 
transducer. Figure 2 presents a data set of K Q versus AP in the 
linear Stokes-Darcy regime at room temperature 0 . 

Additional results are obtained readily using flow measurements of 
Newtonian fluids. Figure 3 shows quasi - steady outflow of liquid 
from an elevated reservoir through a wire-m-tube system /7/ into 
an outer bath. In the laminar range, the level difference decays 
exponentially with time t : Az/Az Q = exp (- t/T Q ); T c character - 
istic time . The latter is proportional to the reciprocal permea - 
bility, to the length, and the cross sectional area ratio, to the 
kinematic viscosity , and to g -1 ; (g gravitational acceleration) . 

Another method is the utilization of "pseudo-classical" flow in the 
liquid breakthrough mode of vapor -liquid phase separation (Fig. 4). 
He II is kept on the downstream side, however high pressure causes 
He I conditions on the upstream side /8/. The resulting Revalues 
versus the externally applied AP are shown m Figure 4. 


DATA DISCUSSION AND COMPARISON WITH OTHER RESULTS 


The temperature influence on the permeability during He^ gas flow 
is displayed m Figure 5 as volumetric flow rate(V) versus AP . Low 
K» ,-values are indicated. There is a transition to non-linearity 
with even lower K gf ^-results . The insert of Figure 5 shows Kj. , low 
values of K -- caused by low-T walls effects on thermometry, and 
high K e £--values for "high" -T walls . 

Compressible flow effects appear to be involved implicitly in the 
data of RM German /9/ for stainless steel plugs. Figure 6 displays 
equivalent Ergun diameters (D-) as a function of porosity. The Darcy 
permeability of Ergun' s equation is 


(K D 5 E 


e 3 /[ (W) 2 iso ] 


(5) 


the bed of 


(D mean particle diameter of 
If*the sintered plug has similar properties, there should be an equi- 


valent Ergun diameter 

V*: 

D_ = L 
r. c 


near-spherical particles) . 
e 

relatively close to 
(1- e )/ (e 3 /lS0) 1/2 (6) 


The D E -values at porosities near closely packed spheres appear to 
vary But little with £ . Thus, m the vicinity of e» 1/3 , D £ ought 
to be a reasonable measure of the sintered particle system. At high 
and low porosities flow patterns are at variance with the near- 
spherical, closely packed particle system. It is noted that com - 
pressibility affects the high K ef ^-aata / / in Figure 6. Upstream 
pressures near 4 bar , with downstream pressures near 1 bar, are 
probably causing local sound speeds. K - r is less likely represent- 
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ative of Darcy convection once Mach number effects are omitted. 

Data are compared readily as long as geometric and flow similarity 
is maintained. A modified diameter D E e /(I -£) is used . Fig. 7 
presents this diameter versus L c . Exact monosphere assemblies 
ought to be described by a unique point in this diagram. The data 
spread, for a specified plug, is in part related to tolerances of 
the manufacturing process. It is seen that various plugs are rea - 
sonably separated from each other. Figure 8 displays the Dg - dis- 
tribution of the present plugs with 2 pm limiting pore size , and 
minimum particle retention respectively (manufacturer quotation ) . 
Two plugs of one manufacturer are lower m D £ than the remainder 
of another manufacturer. 

The possibility of Knudsen transport contributions has been 
checked by measuring the permeability with N 2 gas. For a 2 urn 
plug, K-. has been found to be 10 % larger than the He^ gas value. 
This difference, translated into Dj;, is small compared to the width 
of the Dg-distnbution in Figure 8. 

CONCLUSIONS 

From the present results it is concluded that the Darcy permeabil- 
ity is quite useful for porous media characterization. An equiva- 
lent diameter may be deduced readily, once K D and are known. 
Precautions however are required to avoid disturbing effects. 

Concerning low-T methods , we note that several He II-based tech- 
niques are useful, e.g. transport near the critical Leiden veloc- 
ity, zero net mass flow and the low-T limit of vapor-liquid phase 
separation. In addition, the acoustic index of refraction may be 
potentially useful /9/ . 
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Fig. 5. Volumetric flou rate V vs. 
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Fig. 6. Ergun diameter of equivalent 
near-spherical particle 
assembly; (German / 6 /) ; 
Material :Stainless steel; 

• Darcy permeability , air; 

Q Data for large AP ; 

Insert: Porosity vs. mean 
■ particle diameter/b/. 
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THERMODYNAMIC PERFORMANCE EVALUATION OF LIQUID PUMPS BASED ON THE FOUNTAIN EFFECT 
OF SUPERFLUID LIQUID HELIUM II 


S.W. K. Yuan, T.H.K. Frederking and R.M.Carandang 
Universlcy of California , Los Angeles 90024 CA USA 


Thermodynamic performance of thermally activated liquid pumps 
("thermopumps") for He II has been evaluated for a heater-acti- 
vated system("Leiden pump") and a cooler-activated device(Klip- 
p/ng type pump"). The theraopump wit h a cooler is found to be 
inferior to the Leiden pump from the point of view of refriger- 
ation power requirements, however it is attractive when cheap 
refrigeration is availaole. Transport rates are conveniently 
based on classical Darcy permeabilities (Kjj) when porous plugs 
are used. Examples are power law approximations for the non - 
Newtonian flow describing the ratio of the He II permeabil- 
ity Kjj to K_. Experiments have been conducted with a Leiden 
pump fitted with a stainless steel plug characterized by 
a nominal size of 2 pm . 


INTRODUCTION 

Recenc studies of forced flow of liquid He II , e.g. [l, 2], raise questions con - 
eeming the availability of suitable pumps. Several mechanical pumps are potentially 
useful, e.g. [ 3, 4], also for He II. However in some cases, pumps without moving 
parts may give great reliability. These types of pumps are available in the He II 
temperature range (T<T^ = 2.172 K) on Che basis of fountain effect utilization. The 
pressure difference APf’ , generated by a temperature difference. Is established 
across a porous medium. The latter is to suppress backflow of heat and normal fluid 
of He II toward the low cemperature side of the plug. Simultaneously a reasonable 
throughput of liquid through the pump is to be established. There are two proto - 
types of thermopumps available: The heater-activated pumo is called a Leiden pump 
because of extensive work at Leiden with this type of pump. The cooler-activated 
Pump, so far, appears to be operated in batch modes for liquid vessel replenishment 
by Klipping ec al. [ 5] . It is the purpose of the present scudies to evaluate 
thermodynamic work requirements and Co assess transport conditions. Experiments have 
been conducted with a stainless steel pump of the Leiden type. 


THERMODYNAMIC EQUILIBRIUM STATE CONSTRAINTS 

It is assumed that a near-saturated He II bach is available ac the specified cemper- 
ature T (= initial T-value of the Leiden pumo , or final T-value of Che Klipping 
Pump) . In the limit of a small temperature difference AT « T , the ideal additional 
refrigeration energy required per thermal energy of the pump is 

AS/Q “ T AT / T 2 , AT / T « 1 (1) 

e 

The uork is based on the assumption chat the environmental temperature (T g ) of aoout 
300 K is the ulcimace thermal energy rejection Cemperature. The puraomg conditions 
3re skecched in Figure 1 in the pressure-cemoeracure diagram P(T); (P v vaoor pres - 
Sure curve). The laeal Klipping pump uses a low ulcimace temperature with a small 
Entropy (p s) per unit volume ,( p liquid densicy) . The fountain pressure difference 
is /T+AT 

A? t •/ p S (T) dT (2) 


+) Presented at 10th International Cryog. Engng. Conf. ICEC-10, 
Helsinki, 1984; (paper 0fe-5) . 
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The insert (a) of Figure 1 shows a simple scheme Eor cne Klipoing pump, and insert 
(b) sketches a Leiden pump system with an aftercooler(AC) . In the subsequent dis - 
cussion a specified pressure difference A? is assumed in order to compare pump 
performance quantities. Thus, the pumps in Figure 1 supply the same AP.j.» AP , and the 
ultimate pressure is snifted by a constant P-change with respect to the vapor pres- 
sure curve P y (T). The Klipping pump is characterized by a required AT^, and the 
Leiden pump needs AT^p . As the entropy is a monotonieally increasing tunction of 
T, Equation (1) implies a larger AT for the Klipping pump than for the Leiden 
pump. In the Klipping pump system (insert a), a liquid precooler (PC) is operated 
by means of a refrigerator (KR). A constant mass flow rate (m) is supplied from a 
He II bath (RB) during steady flow of liquid through the porous plug (SP) of the 
pump. The Leiden pump (insert b) has a heater (H) downstream of the plug (SP). In- 
troducing a power law approximation for the entropy, S ~ T 04 , one obtains thermal 
energy needed for each thermopump as Q *<k(AP_ /<f ). Thus, the ideal thermodynamic 
work per mass is obtained by insertion of Q into Eq.(l) for AE/Q. For an arbi- 
trary AT the work per mass of the Klipping pump is obtained as 


W/m = ^ ( AP 7 /<f) [(T - AT) 


-1 


_1 ] T 


(3) 


In the related (W/m)-equation for the Leiden pump, a lower numerical value results 
because of a generally higher temperature. Subsequently, a particular cooler system 
option is considered for the Klipping pump (insert of Fig. 2 ). A vacuum pump at 
an elevated temperature is to exert T-control via vapor pressure control. A mass 
flow rate m^ is taken from the main bath (RB) to an auxiliary bath (AB) at (P V *4P ) 
whose vapor is pumped off by the vacuum pump (VP). The ideal isothermal work of the 


vapor removal system is considered for comparison 


(R gas constant) 


W KP /a KP ■ R T j [la P atm " ln(P v " ^ v >] (4) 

Figure 2 presents the work of^Equation (4) for various T.-levels below T e> i.e- for 
ideal cold compression of He . In the lower part of Fig.l the ideal W/m of Eq. 
(3) is included. For equal flow rates of steady operation, and at T below 4 K, the 
effectiveness of the pump, expressed as W^/W > in the best case does not exceed 
10 Z significantly near 1.9 K. For cold compression above 20 K, the order of 1 Z 
is reached, and at T g the pump effectiveness becomes quite small. 


TRANSPORT CONDITIONS 

The thermopump does not only have to satisfy thermostatic requirements. In addition 
the mass flux in conjunction with the required plug cross section has to be satis- 
factory. One may adopt a simplified classification of flow with the following three 
regimes : 1. Superflow (Landau regime with zero flow resistance);!. Transition re- 
gime with vortex shedding processes leading to vortex flow resistance ; 3. "Classical" 
Newtonian fluid flow. Near-critical operation at vortex shedding onset has a small 
flow resistance. For a tie-in with Darcy convection of Newtonian fluid, the recipro- 
cal resistance ratio is considered a permeability ratio (Kj^/Kp) for the present 
studies ; (Kjj Darcy permeability). In the vortex flow resistance regime the pover 
law approximation of Vote et al. [6] is adopted : For instance, for a sintered 
1 pi plug (nominal size), the flow may be described by 

(Ku'V ‘ ^CPsVK P S i 7 T l> ? L c 3 /> {n 2] " 3/4 (5) 

(L^ = Kj ' , normal fluid shear viscosity , ^ ■ 0.36 + 0.04) . The oermeability 
ratio of Eq. _(5) is rewritten in terms of the volumetric flow race per Jo£al PJ- U S 
cross section (v Q ) making use of the power law [6] with v o « ( p S |7Ti) 

“JV 4 <?.'?>* ' <e\> L c '*w> 3 (6) 

According to the power_law aoproxmation , there is a discrete value of the super- 
fluid transport race (v p> /p g ), multiplied times L^ ' , at the ooservacion limit 
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for AP and AT respectively; (Leiden criterion for the superfluid "critical" veloc- 
ity). At low speed there are departures fron the power law, as K grows beyond all 
Units ac the cessation of the Landau regiae. At high speed the vortex shedding 
processes become disordered and chaotic. Figure^ 3 shows the permeability ratio. 
Equation (6) , versus the normalised flow rate (v Q p L /v^) ■ (? s v os I* / Ij^) at sever- 
al temperatures for the plug of Ref. [6] omitting tfie viciniey°of the lambda 
temperature . In the frame of Equation (6), flow quantities are affected sensitive- 
ly by data uncertainty. Therefore, a variation of the constant has been incorpo- 
rated in the data comparison (caption of Fig. 3). 

EXPERIMENTS 

A small chermopuap of the Leiden type has been used for the purpose of replenishment 
of a liquid jackec simulating a liquid He II vessel during vapor - liquid phase sep- 
aration experiments [7], A 2 pm stainless steel plug with a diameter of 1.27 cm 
and a thickness of 0.3173 cm has been employed with a soldering joint between a sur- 
rounding stainless steel tube. The pump is shown schematically in the insert of 
Figure 4. Satisfactory performance was observed. In line with the experimental pro- 
cedures used, [7], a drift in temperature was permitted (Fig. 4) . In Figure 4 the 
liquid level z is shown as a function of the time t . In addition the speed v Q ~ 
dc/dt is displayed. 

CONCLUSIONS 

Thernopumps have advantages when no moving parts are desired for enhanced reliabil- 
ity.when the space required for the pump is noc a crucial factor, and when the nec- 
essary AP-values are noc excessive. The effectiveness may be low with some pumping 
schemes, and the Klipping type pump, ac a first glance, appears Co be inferior to 
the Leiden pump. This is not the case however for large scale systems with suffi- 
cient refrigeration available at suitable T-levels. 
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Fig. 1. Pressure-temperature diagram with 
^ state changes in thermopumps; 

Insert a: Klipping type pump for 
steady flow of liquid; 
Insert b: Leiden pump with heater 
and after cooler. 
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Fig. 2. Comparison of energy requirements 

for AP - 100 mbar; ^KP^ m KP> 

Insert : Example of Klipping 
type system 
(schematically) . 
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Fig. 3. Permeability ratio for various 

bath temperatures [6] of a 1 pm 
plug with Kj* 2.5-10 cm2 . 


. Data of tnermopumps experiments 
versus time ; KL=3.4 10~® cm2 
at room temperature, 

Insert : Pump system ( H heater, 
SP porous plug,TH thermometer; 
(z liquid level position in 
vessel) . 












This report outlines work performed during Summer 1984 
(from June to Sept. 30) in the area of porous media for use in 
low temperature applications. Recent applications are in the area 
of vapor - liquid phase separation, pumping based on the fountain 
effect and related subsystems . Areas of potential applications 
of the latter are outlined in supplementary work. 

Experimental data have been taken to the extent permitted 
by the limited funding . The linear equations of the two-fluid 
model are inspected critically in the light of forced convection 
evidence reported recently. It is emphasized that the Darcy 
permeability is a unique throughput quantity in the porous media 
application areas whose use will enable meaningful comparisons 
of data not only in one lab but also within a group of labs doing 


porous plug studies. 


